##funcao 
 
LP <- function(data) {
  N <- c()
  for (i in 1:length(unique(data$inq))) {
    tmp <- subset(data, inq == i)
    N[i] <- with(tmp, ( (sum(n1) + 1)*(sum(n2) + 1) / (sum(m2) + 1) ) - 1)
  }
  x11()  # abre um dispositivo grafico
  par(mfrow = c(3,1))
  plot((N),main=("Plot de Lincoln Peterson"), xlab="inquerito")
  barplot((N),main=("Barplot de Lincoln Peterson"))
  boxplot((N),main=("Boxplot de Lincoln Peterson"))
  return(N)
}


##Testando a funcao com example 


####Banco de dados creado inventado para simular os dados de uma populacao x 

tabela <- matrix("Nas",nrow=200,ncol=4)
idade <- sample(c("j","a","i"),prob=c(0.25,0.50,0.25),size=200,replace=T)
sexo <- sample (c("m","f"), prob=c(0.60,0.40),size=200,replace=T)
condrep <- sample (c("c","e"),prob=c(0.30,0.70),size=200,replace=T)

for(i in 1:200)
{
  tabela[i,1] <- i
  tabela[i,2]<- idade[i]
  tabela[i,3]<-sexo[i]
  tabela[i,4]<- condrep[i]
}
tabela <- data.frame(id=tabela[,1],idade=tabela[,2],sexo=tabela[,3],con_rep=tabela[,4])
tabela


####Creacao de um inquerito simulado composto por 20 eventos de captura e recaptura 
####como entrada para provar a funcao

inquerito <- NULL
for (i in 1:20) {
  ran <- sample(200, 1)
  tmp <- sample(tabela[, 1], ran)
  tmp <- tabela[tmp, ]
  inq <- cbind(tmp, sample(c(0, 1), ran, replace = T),
               sample(c(0,1), ran,replace=T),
               sample(c(0,1), ran,replace=T),
               rep(i, ran))
  
  inquerito <- rbind(inquerito, inq)
}
names(inquerito) <- c('id', 'idade', 'sexo', 'con_rep', 'n1', 'n2', 'm2', 'inq')

inquerito

###Testando!!!

LP(inquerito)
  

 

